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Abstract. Measurements of a variety of hadron species in pp and Au+Au colUsions 
at 200 GeV are presented and studied within the framework of chemical and local 
kinetic equilibrium models. The extracted chemical and final kinetic freeze-out 
temperatures and collective flow velocities are discussed as function of centrality. The 
results suggest that Au+Au collisions of various centralities at RHIC always evolve 
toward the same temperature at chemical freeze-out, followed by cooling and expansion 
toward kinetic freeze-out. 



1. Introduction 

It is predicted that at sufficiently 
high energy densities matter is no longer 
confined to hadrons but exists in a 
deconfined state of quarks and gluons, 
the quark-gluon plasma (QGP). This 
new state of matter is expected to 
form in relativistic heavy ion collisions 
at RHIC. QGP signals may remain in 
the bulk properties of the collision at 
final state. For example, the final 
collective transverse radial fiow, due to its 
cumulative nature, contains any collective 
fiow generated at partonic stages, and 
particle yield ratios may refiect statistical 
nature of the hadronization process. We 
present here measurements of a variety 
of hadron species (vr^, K^, K^, K* ^ 
0, p, p, A, A, S, H, Q + Q) in pp 
and Au+Au collisions at 200 GeV by 
the STAR experiment, and a detailed 
investigation of the final hadronic state pi 
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Figure 1. 

Identified particle spectra in 5% most central 
Au+Au events measured by the STAR experiment. 
Black solid lines show blast-wave fit results for 
charged pions, kaons and protons. Dashed 
lines represent blast-wave prediction for multi- 
strange baryon spectra shapes based on freeze-out 
parameters of tt, if, p. Red lines are blast-wave fits 
for S and Q, spectra. 

perties of such collisions in the framework of 
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chemical and local kinetic equilibrium models. 
2. Data Analysis 

The STAR detector is described in [T]. The main detector, Time Projection Chamber 
(TPC), tracks charged particles in a uniform magnetic field of 0.5T. vr^, K^, p, and p 
are identified be ionization energy loss; K^, K*, 0, A, A, S, H, and f2 + Q are identified 
by their decay topology. Corrections for tracking inefficiency, detector acceptance as 
well as decay losses were applied. Pion spectra were corrected for weak decay feed- 
down, muon contamination and background pions; all other spectra include feed-down 
from weak decays. Additional corrections for p+p spectra were applied to account 
for inefficiency and false rate of the vertex finder. Minimum bias Au+Au events were 
subdivided into centrality classes based on measured number of charged particles within 
the pseudo-rapidity range |?7| < 0.5. 

Figure 1 shows comparison of the 
spectra for different particle types in cen- 
tral Au+Au collisions. Heavier parti- 
cles show strong ffattening of the spec- 
tra, consistent with radial ffow resulting 
from high pressure gradient in the reac- 
tion zone. We fit n"^, K^^ p, and p 
spectra with a hydrodynamically moti- 
vated model with common flow velocity 
(also known as 'Blast-wave' model) with 
kinetic freeze-out temperature T^in and 
mean flow velocity (/?) as flt parameters 
[21 . Fit results are shown in Figure 1, 
demonstrating good description of vr, K, p 
spectra by such a fit (solid lines). Dashed 
lines shown in Figure 1 illustrate that blast-wave model fit to the common particle 
spectra fails to reproduce spectra shape of rare particles. 

Blast- wave fits for rare particle spectra are shown as red lines in Figure 1. Extracted 
values of Tkm and are given in Table 1. It has been suggested |3] that multi-strange 
baryons freeze-out earlier than the common particles, possibly due to their smaller 
interaction cross-section with the bulk of the collision zone. 

Integrated dN/dy values are extracted from the fit. Ratios obtained from those 
yields are fitted within the framework of statistical model with 4 parameters: chemical 
freeze-out temperature (Tc^), baryon and strangeness chemical potentials (/i^, /is), and 
strangeness suppression factor (7^) |4j|. Points shown in Fig. 2 represent measured STAR 
data for the most central events, blue lines are statistical model results based on fit to 
the vr, p measurements. The obtained values of fit parameters are listed in Table 
2 (first column). As seen in the plot the relative abundances of strange particles are 
underestimated. Results of statistical model fit to all the ratios (data available are listed 
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Figure 2. 

Measured particle ratios (symbols) and statistical 
model predictions (lines). Blue lines represent fit 
results based on tt, K and p data. Red lines show 
results based on tt, K, p, A, S, fl and (j> data. 
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in the Table 3) are shown in red in Fig. 2. The thus obtained values of Tch, /is, and /i^ 
agree with the vr, K, p fit results, while strangeness suppression factor is larger, giving 
better description of strangeness production. The results are given in the last column of 
Table 2. On the other hand both fits fail to reproduce STAR resonance measurements, 
which may be due to rescattering/regeneration effects after chemical freeze-out. 

Open symbols in Fig. 3 summarize the centrality dependence of freeze-out 
parameters extracted from tt, p measurements: Tch ~ 160 MeV and is independent 
of centrality, Tkm decreases from ^ 140 MeV (peripheral) to 90 MeV (central) and 
(P) and 7s increase with centrality 0. The results suggest that Au+Au collisions of 
different initial conditions always evolve to the same chemical freeze-out condition, and 
then cool down further to a kinetic freeze-out dependent of centrality. The expansion 
of the system gives rise to collective flow. Significant drop in the temperature suggests 
relatively long time duration between two freeze-outs, estimated to be at least 6 fm/c 
for central events [H]. 

It also appears that the chemical 
freeze-out temperature from our fit is 
close to the Lattice QCD predicted ^ 
critical temperature of ~ 170 MeV 
(shown red line in Fig. 3). 

Figure 3 also shows the kinetic freeze- 
out parameters from blast-wave fits to 
other particles. The S results of those 
fits are shown as function of centrality 
(red symbols); the fit results of other 
particles are for central collisions. The 
data, although with large uncertainties, 
seem to indicate sequential kinetic freeze- 
out of particles: and S decouple from 
the system right after chemical freeze-out, 
followed by other particles (A, K*, tt, 
p). Also it appears that the kinetic freeze- 
out of S" and S"*" baryons coincide with 
chemical freeze-out for all centralities. 
It is therefore a tempting to conclude 
that the {(3) = 0.45 ± 0.1 c extracted 
from H~ and spectra characterize 
radial flow at chemical freeze-out, which 
is possibly build up prior hadronization in 
the earlier, probably partonic, stage of the collision. 
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Figure 3. 

Upper panel: chemical and kinetic freeze-out 
temperatures from tt, K, p data; kinetic freeze-out 
temperatures for S, fl, A and K* data separately. 
Lower panel: Mean transverse flow velocity for tt, 
K, p (black), S (red), (green), A (blue), and K* 
(purple). The fit results for pp data are shown at 
dNch/dT] ~ 0. 
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3. Conclusions 

In summary, STAR has measured a wide variety of hadron species produced in 200 GeV 
Au+Au collisions. The data are studied within the framework of chemical and local 
kinetic equilibrium models to investigate final hadronic state properties of heavy ion 
collisions at RHIC. The major results of our study are: 

• A chemical equilibrium model fit to the ratios yields a Tch insensitive to centrality. 
The estimated value of Tch ~ 160 ± 6 MeV is close to the QCD predicted phase 
transition temperature. 

• Kinetic freeze-out temperature of multi-strange baryons shows no sensitivity to 
collision centrality and coincides with chemical freeze-out temperature for all 
centralities. The obtained flow velocity {(3) = 0.45 ± 0.1 c may be considered 
as that at hadronization. 

• Blast wave model fit to the transverse mass spectra of different particle species 
indicates likely sequential kinetic freeze-out (51, A, 7T,K,p,K*) and allows to 
estimate radial flow velocity at the chemical as well as kinetic freeze-out. 

• Kinetic freeze-out parameters obtained from blast-wave fits of tt,K^p spectra are 
correlated with centrality: the more central the collision, the lower extracted 
temperature and the higher collective radial fiow velocity. For the 5% most central 
collisions, Tkin = 89 ± 10 MeV, (/?) = 0.59 ± 0.05 c. The drop in temperature 
from Tch to Tkin and the development of strong radial fiow suggest a significant 
expansion and long duration from chemical to kinetic freeze-out, estimated to be 
at least 6 fm/c in central collisions. 

The results seem to suggest the following picture at RHIC: Au+Au collisions of 
various centralities, despite different initial conditions, always evolve toward the same 
chemical freeze-out temperature (or possibly hadronization); chemical freeze-out is 
followed by cooling and expansion, and a sequential decoupling of particles dictated 
by their hadronic cross-sections: multi-strange baryons decouple at chemical freeze-out, 
followed by other particles decoupling at later times. 
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Table 1. Blast- wave model fit results for 5% most central Au+Au collision. Flow 
profile parameter n is free parameter for tt, K, p fit. The resulting value of n = 0.82 
is fixed for rare particle fits. 



Particle 



TT, K, p 

K* 
A, A 



kin 



(MeV) iP) (c) 



89 ±10 
75 ±35 
115 ±20 
161 ±20 
179 ± 60 



0.59 ±0.05 
0.62 ±0.05 
0.54 ±0.05 
0.46 ±0.10 
0.45 ±0.10 



Table 2. Chemical frcczc-out paramotors from statistical model fits to the particle 

ratios measured by STAR in central Au+Au collision. 



Parameters t^,K ,p ratios fit i: ,K ,p,K,(t>,^,Q, ratios fit 



Tch (MeV) 
MB (MeV) 
Ms (MeV) 

Is 



157 ±6 
22 ±4 
3.8 ±2.6 
0.86 ±0.11 



160 ±5 
24 ±4 
1.4 ± 1.6 
0.99 ±0.07 



Table 3. Particle ratios (errors are systematic) and statistical model fit results for 5% 
most central Au+Au collision. Ratios from the bottom portion of the table were not 
included in fit. 



Ratio 


STAR data 


Model prediction 




1.01 ±0.02 


1.00 


K-/K+ 


0.96 ± 0.03 


0.93 


p/p 


0.77 ±0.04 


0.77 


K-J-K- 


0.15 ±0.02 


0.17 


p/n- 


0.082 ± 0.009 


0.091 


A/tt- 


0.054 ±0.006 


0.052 


A/tt- 


0.041 ±0.005 


0.043 


S-/7r- 


7.8 -10-3 ±10-3 


6.0-10-3 


S+/7r- 


6.3 • 10-3 ± 8 • 10-^ 


5.5-10-3 




9.5 • 10-"* ± 10-4 


10.1 • 10-^ 


fi/n 


1.01 ±0.08 


0.95 


A(1232)/p 


0.24 ±0.04 


0.12 




0.12 ±0.03 


0.14 


K*JK- 


0.205 ± 0.033 


0.329 


A(1520)/A 


0.033 ±0.017 


0.061 



